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SHOCK  CURVATURE  AND  GRADIENTS  AT  THE  TIP  OF 
POINTED  AXISYMMETRIC  BODIES  IN  NONEQUILIBRIUM  FLOW 

ABSTRACT 

The  shock  curvature  and  flow  variable  gradients  at  the  tip  of  a 
pointed  body  caused  by  nonequilibrium  effects  are  considered. 

Coordinates  introduced  by  Chester  are  used  since  they  offer  a  convenient 
way  of  treating  the  boundary  conditions.  The  desired  functions  are 
obtained  by  solving  numerically  a  system  of  linear  ordinary  differential 
equations.  These  equations  have  a  singularity;  the  nature  of  the  singu¬ 
larity  is  found  analytically,  and  its  numerical  treatment  is  discussed. 
The  specific  nonequilibrium  effect  considered  is  vibrational  relaxation 
in  a  pure  diatomic  gas.  Representative  results  are  given  for  flow  of 
Ng  over  a  cone  for  a  comprehensive  range  of  Mach  number  and  cone  angle. 
There  is  a  point  analogous  to  the  Crocco  point.  These  results  compare 
favorably  with  those  obtained  by  South  and  Newman  using  an  approximate 
method.  Another  check  is  made  by  comparison  with  characteristic 
calculations  extrapolated  to  the  origin. 
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—  angle  between  x  jaxis  and  tangent  to  shock  wave 

/ 

—  c  /c  v  =  1.4  in  present  calculation 

Pa  a 

_  0  for  ideal  gas  flow,  =  1  for  nonequilib.  flow 

—  Chester  variable  £  =  2 i[i'/[p  'q  'y/2J 

CO  CO 

—  Chester  variable  Tj  =  y 

—  angle  between  velocity  vector  and  x'  axis 

—  characteristic  vibrational  temperature  of  the  gas 

—  density  p  =  p '/p  ' 

— •  vclflYfl+.inri  f.nmo  T  ~  t^/t  1 
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'  b 

—  stream  function 

—  function  defined  in  equation  (2.6) 

—  OE/ari)^  =  0 
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temperature 
(du/ftT))^  =  Q 

(^v/ciTl)^  =  Q 

0  Vt  ' 
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LIST  OF  SYM30IS  (Contd) 


Subscripts 

oo  —  free  stream  conditions 

b  —  body 

w  —  shock  wave 

h  —  refers  to  homogeneous  differential  equations  and  their 

solution 

n  —  refers  to  particular  solution  to  nonhomogeneous 

differential  equations 

Superscript 

—  denotes  quantity  with  dimensions 


1.  INTRODUCTION 


When  nonequi librium  effects  are  considered  in  the  fluid  dynamic 
equations,  the  classical  similar'  solutions  describing  supersonic  flow 
over  wedges  and  cones  no  longer  exist.  To  gain  s’ome  insight  into 
nonequilibrium  effects  one  can  seek  solutions  valid  in  restricted  regions 


x  "L  . 


of  a  flow,  e.g.,  the  vicinity  or  tne  tip  of  a  wedge  or  cone.  The  first 
approximation  to  the  flow  here  is  the  frozen  (similar)  solution.  The 
next  approximation  brings  into  evidence  the  nonequilibrium  effects,  such 
as  the  curvature  of  the  shock  wave.'  The  purpose'  of  this  work  is  to  find' 
this  approximation  for  the  case  of  flow  over  a  cone  and  pointed  ogives 
in  general. 


i _ . 


ine  corresponding  jjruuiem  iOr 


1  t  j-m  s— . 

wcagc 


^  ucmicj  , 


i  r\C  i  ’i  A  , 


by  explicit  solution  of  algebraic  equations.  For  cone  flow,  differential 
equations  with  two-point  boundary  conditions  (at  shock  and  body  surface) 
must  be  ^solved.  The  coefficients  of  these  equations,  which  involve  the 
frozen  conical  flow  solution,  give  rise  to  a  singular  point,  at  least 
in  all  choices  of  coordinates  investigated. 


There  is  a  close  analogy  between  the  problem  of  approximating  the 
nonequilibrium  flow  in  the  neighborhood  .of  the  tip  of  a  wedge  or  cone 
and  that  of  approximating  the  equilibrium  or  frozen  flow  near  the  tip 
of  a  pointed-ogive  body.  The  latter  problem  was  solved  by  Crocco  (1937) 
for  two  dimensional  flow  and  by  Shen  and  Lin  (1951)  for  axisymmetric 
flow.  In  order  to  treat  the  singularity  for  the  axisymmetric  case,  and 
especially  its  effects  on  higher  approximations,  Kogan  (1956)  formulated 
the  problem  in  terms  of  Crocco’s  stream  function.  This  refinement, 
however,  is  not  included  in  the  values  of  shock  curvature  hitherto 
published,  Shen  and  Lin  (1951)  ,  and  Bianco  et  al,  (i960). 

In  the  method  of  solution  adopted  here  for  determining  shock 
curvature  and  gradients  caused  by  nonequilibrium  effects,  these  same 
quantities  for  frozen  flow  over  ogival  bodies  are  obtained  as  a 
by-product.  The  specific  source  of  the  nonequilibrium  effects  considered 
here  is  vibrational  relaxation  of  a  diatomic  gas;  however,  the  present 
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method  can  be  employed  for  dissociation,  ionization,  etc.  The  model  of 
a  vibrationally  excited  gas  serves  as  well  as  any  other  example  to 
illustrate  the  technique  and  effects. 

A  distinctive  feature  of  this  analysis  is  the  choice  of  coordinates. 

Specifically  the  coordinates  used  are  those  introduced  by  Chester  (1956) 

in  his  study  of  hypersonic  flow  over  blunt  bodies.  The  chief  advantage 

of  this  choice  is  the  ease  of  handling  the  boundary  conditions  on  the 

shock;  for  here  the  unknown  shock  curve  is  mapped  into  one  of  the 

coordinate  lines.  Also  the  singular  point  at  the  tip  is  spread  out  into 

i 

a  line  interval  along  the  second  coordinate  axis.  The  latter,  but  not 
the  former,  advantage  is  also  gained  with  polar  coordinates.  The 
differential  equations  now  have  an  integrable,  inverse  one-half  power- 
type  of  singularity  at  the  origin.  This  could  be  handled  in  the 
standard  way  —  using  expansions  to  match  the  numerical  solution;  instead 
the  entire  problem  is  treated  numerically,  the  authors  believe  this 
procedure  to  be  more  efficient.  The  singularity  is  treated  by  a  "detailed 
approach  to  the  body"  which  involves  decreasing  the  step  size  in  the 
Runge-Kutta  numerical  integration  procedure. 

The  motivation  for  this  work  arose  in  the  computation  of 
nonequilibrium  flows  over  cones,  Sedney  and  Gerber  (1963),  by  the  method 
of  characteristics.  One  is  forced  to  initiate  that  calculation  by 
assuming  s.  fin  its  frozen  flow  r,0gi  on :  Tiic  of  feet  of  "this  error*  esn  i>e 
relatively  large  and  propagate  two  or  three  lengths  of  the  original 
frozen  region  downstream  even  though  the  grid  size  is  twenty  to  thirty 
times  less  than  this  length.  The  behavior  near  the  tip  can  be  studied 
by  the  laborious  process  of  choosing  successively  smaller  frozen  flow 
regions . 

The  results  of  the  present  paper  can  he  used  as  a  check  on  those 
from  the  characteristic  method.  Alternatively,  they  can  be  used  as 
input  to  start  the  characteristic  calculation.  With  the  improved 
approximation  to  the  flow  near  the  tip,  a  larger  grid  size  is  possible 
and  a  saving  on  computation  time  realized. 
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2 .  FLOW  EQUATIONS  ' 

The  axi symmetric  flows  to  be  ’considered  are  steady,  inviscid,  and 
« 

isoenergetic.  The  equations  for  conservation  of  mass ,  momentum,  and 
energy  are ,  in  non-dimensional  form. 


fe  (ypu)  +  (ypv)  =  0 


(2.1) 


yM  2  fu  +  v  &]  =  -  - 
“  L  3yJ  p  3x 


(2.2) 


„2  r  3v  9vl  1  3p 

”  L  3x  3yJ  p  3y 


(2.3) 


^  M  2  (u2  +  v2)  +  T  +  E  =  h,  (2.U) 

d  °°  "C 

See  Figure  1  for  notation  (where  primes  indicate  dimensional  quantities); 
p  is  density,  p  is  pressure,  T  is  temperature,  E  is  vibrational  energy 
per  unit  mass,  y  is  the  ratio  of  frozen  specific  heats,  and  h^  is  the 
constant  total  enthalpy.  Variables  are  made  dimensionless  as  follows: 

lit 

lengths  by  q^  where  r  is  the  relaxation  time  evaluated  at  the 

frozen  conditions  on  the  body  tip  and  q^  is  free  stream  velocity; 

velocities,  pressure,  density,  and  temperature  by  their  free  stream 

/  / 

values;  vibrational  energy  by  c__  T  ,  the  frozen  enthalpy  in  the  free 

pa  CO 

stream.  The  model  chosen  for  the  diatomic  gas  with  vibrational  relaxation 
is  discussed  in  Sedney  (1961).  The  perfect  gas  law  is  assumed 


p  =  pT 


(2.5) 


and  the  rate  equation  is 

u  H  +  v  -^  =  [E*  (T)  -  E]/t  =  $ 


(2.6) 
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Cross-sectional  view  of  flow  field  in  the  physical  pi 
(qf  =  flow  speed;  u',  v'  =  velocity  components  in 
coordi nate  di recti ons .)  i 


t 


SHOCK 


where  x  is  the  relaxation  time  referred  to  x  y  ,  and  E*(T)  is  the  local 
equilibrium  value  of  E,  given  by 

E*(T.)  ='(2Z/7)/(eZ/T.-  1) 


with  Z  =  f?  '/T  '  ,  £5  •  =  characteristic  vibrational  temperature.  The 

:  v  00  v  .  - .  - 

variation  of  x  with  T  and  p  is  approximated  by 


t  =  (p^/p)  exp  [(c/^)1/,3  -  (c/T^)1^3] 


and  the  subscript  b  refers  to  conditions  on  the  body  at  the  tip.  The 
constant  c  is  evaluated  for  various"  temperature  ranges  by  fitting 
experimental  data;  e.g.  ,  those  of  Blackman  (1956)  or  Millikan  and  White 
(1963).  "  '  ' 

The  independent  variables  x,y  are  transformed  to  c,  ti  by 


?  =  2^(x,y)/y  ,  n  =  y  ' 


(2.7) 


where  ^  is  the  dimensionless  stream  function 


3(|j/3x  =  -  ypv,  3tf>/3y  =  ypu 


Thus 


L.=  _2 £iLandL=L+  2.(^_- 


3x  Ti  3?“‘“3y  3n  ' .  n  3  <; 

(2,2).  (2, 3),  and  (2.6)  become,  respectively. 


;  and  Equations'  ( 2 1 )  , 


2  +  i-'/VW 

3?  \vj  3?  \PyJ  3n  \pv J 


(2.8) 


^  =  n  3k’ .  2c  lH 


3n...  ■ 


.  (2.9) 
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SEOHIuOh,  lieeaet 

STTip-  XL 


I 


/<Sr  s_ 

1  2  ]L"  31 


(u2  +  V2)  -  2C  J-  (u2  +  V2 


=  (2-10) 


("If  -2cf)=5-[E*(T)  -E] 


Note  that  if  vibrational  energy  in  the  free  stream  is  negligible 
(E_  =  0) 


ht  =  1  +  (1/2) (y  -  1)  Ma 


(2.11) 


This  would  be  the  case  if  <<  ©v  and  equilibrium  existed  in  the  free 
stream.  If  the  free  stream  is  out  of  equilibrium  E^  will,  in  general, 
not  be  negligible.  This  would  occur  if  a  cone  were  placed  in  the  test 
section  of  a  shock  tunnel  since  the  free  stream  is  generally  frozen  in 
such  a  facility.  In  any  case,  behind  the  shock,  E  takes  its  free  stream 
value . 

The  variables  x,  and  n  were  employed  by  Chester  in  the  study  of 
hypersonic  flow  past  blunt  bodies.  In  the  present  problem  the  shock 
wave  is  attached  at  the  pointed  tip  of  the  body,  but  the  flow  region  is 
mapped  into  the  same  strip  of  the  c,  n  plane  as  for  the  case  of  a 
detached  shock.  This  strip  is  shown  in  Figure  2  where  it  is  seen  that 
the  tip  transforms  into  the  line  segment  n  =  0,  0  <_  £  <_  1.  This 
stretching  is  necessary  to  examine  the  singularity  at  the  tip.  It  is 
also  accomplished  by  using  polar  coordinates  (when  properly  plotted  on 
rectangular  axes);  however,  the  shock  curve  is  unknown  then.  The  same 
strip  is  obtained  if  a  body-oriented  coordinate  system  is  used  and  the 
normal  coordinate  is  normalized  by  the  total  distance  between  the  body 
and  shock. 

If  the  body  is  a  cone  and  the  flow  is  either  frozen  or  in 
equilibrium  so  that  the  flow  is  conical,  then  the  flow  variables  are 
constant  along  the  lines  x,  -  constant.  This  is  shown  using  the  fact 


that  ij)  is  homogeneous  of  degree  two  in  the  variables  x  and  y,  a  result 
which  follows  from  the  definition  of  i|i  and  the  fact  that  pu  and  pv  are 

;l  ’  , 

homogeneous  of  degree  zero.  -The  streamlines  are  always. the  hyperbolae, 
=  constant. 

•*  . 

On  physical  grounds  -one  would- expect  that  the  flow  at  the  tip, 
n  =  0,  should  be  frozen;  this’ is  easily  proved  from  Equations  (2.8)  - 
(2.11).  ..Setting  n  =  0,  all  terms  containing  3/3ri  vanish,  3E/3^  =  0, 
and  the  remaining  equations  give  .  , 


av  =  _  1  [~;bW  ~  g)]/D 

2  l  p3  y 


(2.12) 


(2.13) 


Since  3E/3c;  =  0  along  n  =  0,  E  is  then  constant- there ;  i.e.,  the  flow 

r  -1 

is  frozen  at  the  tip.  The  constant  will  be  the  free  stream  value  of  E. 

The  Equations  (2.12)  are  those  for  conical  frozen  (Taylor-Maccoll) 
flow  in  the  coordinates  n ;  their  numerical  solution  is  a  necessary 
step  in  the  determination  of  the  gradients  since  the  frozen  flow  variables 
appear  as  coefficients  in  the  gradient  equations.  This"  solution  is 

.  a.' 1  ■  '•  . .  ■  •  r  "» 
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obtained  by  integrating  from  the  shock,  t,  -  1,  to  the  body,  5  =  0,  given 
M  ,  y,  and  8,  where  tan  8  is  the  shock  wave  slope.  The  initial  condi- 
tious  at  the  shock  are  the  standard  frozen  shock  relations; 


u  =  1  - 
w 


2  (M  2  sin2  8-1) 
00 

(y  +  1)M^2 


(2. it) 


and  so  forth.  The  half-angle  of  the  body  tip  is  determined  from 

I 

tan  9b  =  (v/u)  =  Q  (2.15) 

One  can  integrate  several  cases  and  then  interpolate  to  find  a  desired 
half-angle  or  set  up  an  iteration  procedure  to  yield  this  result. 


3.  GkADIENT  EQUATIONS 

The  equations  for  the  gradients  of  the  flow  variables  are  obtained 
by  differentiating  Equations  (2.8)  -  (2.1l)  with  respect  to  n,  then 
setting  n  equal  to  zero.  This  requires  the  assumption*  that 
32/3?3n  =  32/3n3C. 

The  following  notation  is  introduced: 


U  =  (3u/3n)n  =  0,  v  =  (3v/3n)n  _  0 ,  0  =  (3E/3n)n  =  Q 

R  =  (3p/3n)  _  n,  P  =  (3p/3n)  _  n 

n  —  u  n  —  u 


(3.1) 


The  set  of  mathematical  conditions  insuring  the  interchangeability  of 
second  derivatives  is  different  from  that  insuring  the  validity  of  an 

2 

expansion  of  the  form  f(r\3z)fo  +  0//3 n^n  +  0( n  Neither  set  can 
be  verified  a  priori . 
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Hereafter,  it  is  understood  that  all  quantities  are  evaluated 
at  n  =  0,  The  differential  equations  for  the  gradient  functions  are 
(where  '  =  d/dc) 


fl  =  (0  -  */v)/(2c)  (3.2) 

U  =  (l/D)  [vG/p  +  h  (C  -  efl)/?]  (3.3a) 

V  =  (1/D)  [-gG/p  +  f  (C  -  efl)/?]  (3.3b) 

R  =  -  [M  2  v?B  +  (g2  +  v2)(A  -  eF)/v  +  j(C  -  e£2)/c1/D  (3.3c) 

00 


where  the  denominator  D  is  the  same  as  in  Equation  (2.13).  The  ■ 
coefficients  g,  b,  f,  and  j  are  known  functions  of  £;: 


g  =  u  -  c/p 

b  =  [pg/.v  -  (y  -  1)  M^2  v?]/(2y  M^p) 
f  =  tc  M^2  (yg  -  u)  +  p]/(2y  M^p) 

■i  =  P  [(y  -  1)  v2  +  g  (y  g  -  u)]/(2y  v)  , 


(3.10 


where  u,  v,  p,  and  p  have  been  obtained  from  the  solution  to  Equations 
(2.12);  •  $  is  given  in  Equation  (2.6).  The  symbol  e  is  equal  to  zero 
for  frozen  flow  (no  vibrational  relaxation),  and  one  for  nonequilibrium 
flow. 

The  coefficients  A,  B,  C,  F,  and  G  are  linear  functions  of  U,  V,  R, 
and  £2: 


A  =  -  M  [y/2  +  (y  -  l)(.up  +  pu)  ]  U  -  (y  -  l)  M  (vp  +  pv)  V 


+  (y  -  i)(pp/p2)R 


(3.5) 
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B  =  (v/v)  U  +  [u/v  +  (cp)/(p2v)  -  2vg/v2]  V 
+  (c/p2)(v/v  +  2p/p )R 

C  =  “  2y£u  -  YC/p)  U  +  Mm2  (v  -  2y£v)  V 

+  (l/p2)(p  +  2 5  p)R  (3.5) 

P  =  p  ^  (l/c)(p/2)(fl  -  */v) 

G  =  pB  +  cA/p  -  e  Upn/p  +  (l/2)(fi  -  $/v)  1 

Equation  (3.2)  can  be  solved  separately,  using  the  initial  condition, 
fiU  =  1)  =  0  (since  E  does  not  change  along  £  =  1) ;  then  a  can  be  treated 
as  a  known  function  of  £  in  Equations  (3.3).  Equations  (3.3)  are  linear 
equations  in  U,  V,  and  R,  homogeneous  when  e  =  0  and  non-homogene ous  when 
e  =  1. 

The  pressure  gradient  function  is  obtained  by  differentiating 
Equations  (2.10)  with  respect  to  n  and  taking  n  =  £  =  0: 

Pb  =  "  Pb  (uU  +  vV)b  (3.6) 

Initial  conditions  are  applied  at  the  shocK  wave  and  £  =  1;  condi¬ 
tions  there  will  be  designated'  by  the  subscript  w.  For  any  function  f 

3f/3n  =  (esc  6)  df/do  =  (esc  B)(df/dg)  K 

=  1 


where  a  is  arc  length  along  the  shock  wave  and  =  dg/do  is  the 
curvature  of  the  shock.  The  initial  conditions  are 


=  (du/de)  esc  8,  =  (dv/dg)  esc  B 

R  =  (dp/dB)  K  esc  8, 
w  w 


(3.7) 
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where  du/d8,  dv/dg5  and  dp/dg  are  obtained  by  differentiating  the 
shock  relations;  e,g.,  Equation  (2. lb). 

The  solution  to  Equations  (3-3)  can  be  written 

U  =  K  11  +  U  ,  V  =  K  V  +  V  j  .R  =  K  R  +  R  (3-8)  - 
wh  n  wh  n5  w  h  n 

where  U^,  V^,  and  R^  are  the  solutions  to  the  homogeneous  equations, 

satisfying  the  initial  conditions,  Equation  (3-T),  with  K  =  1;  U  ,  V 

,i  -  .  .w  n  .  n  ’• 

and  R^  are  particular  solutions  to  the  non-homogeneous  equations  satisfy¬ 
ing  the  initial  conditions  (U  )  =  {V  )  =  (R  )  =0. 

6  n ,  w  n  w  n  w 

With  the  solutions  U,  ,  U  .....  R  determined  it  remains  to  find  K 

h  •  n  n  .  .  w 

for  the  complete  solution.  Equation  (3=8)=  This  is  found  by  applying  the 

terminal  condition.  Equation  (2.15)  to  the  gradient  functions.  From 
Equation  (2.15)  one  obtains 

Vb  =  «b  tan  9b  +  '‘b  (9VSn)n  =  0  .  (3-'9) 


also 


•  (9eb/3n)n  =  o  =  **  csc 


(3.10) 


where  s  d0^/ds  denotes  the'  curvature  of  the  body  at  the 'tip-,  and  s' 
the  arc  length  along  the  body.  Combining  Equations  (3.8),  (3=9)  and 
(3=10)  gives  the  shock  curvature  in  terms  of  the  body  curvature 


K  = 
w 


r  2  -i 

rV  -  U  tan  Gq 

u  sec  0  esc  0 

n  n 

Vu  -  U  tan  0 
L  h  n 

Vv  -  U  tan  6 

L  h  h  ~1 

(3=11) 


Finally,  the  gradients  on  the  body  are  given  by 


(dp/ds)^  =  Pv  sin  0V 


(dE/ds)^  =  Si^  sin  0V 


(3.12) 
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2 


(dT/ds)^  -  (pP  -  pR)b(sin  8b)/Pb 
(dq/ds)^  =  (uU  +  vV)b(sin  6b)/qb 


(3-12) 


where  P  is  given  by  Equation  (3.6). 


Equation  (3.1l)  shows  that  the  shock  wave  curvature  at  the  tip  of 
a  pointed  body  of  revolution  is  the  sum  of  two  terms;  i)  the  curvature 
for  a  curved  body  in  frozen  flow,  and  ii)  the  curvature  for  nonequilibrium 
flow  over  a  cone . 


Previous  calculations,  Shen  and  Lin  (1951),  and  Bianco  et  al.  (i960), 
have  been  made  of  K^/K^  for  ideal  gas  flows  ( e  =  0  in  Equations  (3.3)). 
These  and  calculations  from  the  present  work  indicate  that  K^/K^  becomes 
zero  for  every  axi  symmetric  body  at  some  M  which  produces  a  partially 


subsonic,  partially  supersonic  conical  flow  at  the  tip  (n  =  0).  This 
implies  that  for  each  6^,  there  is  an  for  which,  according  to 
Equation  (3.1l),  I 


(V, 

n 


U,  tan  9), 
n  o 


0. 


This  is  the  "Crocco  point"  in  axisymmetric  flow  (see  Ferri  (1954) ).  It 
is  also  seen  from  Equation  (3.1l),  that  for  these  same  values  of  and 
9^  the  curvature  of  the  shock  wave  in  nonequilibrium  flow  becomes 
infinite  at  the  tip  of  the  cone  (where  1^=0). 

4.  SOLUTION  IN  THE  VICINITY  OF  THE  SINGULAR  POINT 

The  equations  for  the  gradient  functions  (3.2)  and  (3-3)  have  a 
regular  singular  point  at  £  =  0.  Following  Ince  (1926)  the  first  step 
in  examining  the  nature  of  the  solutions  in  the  neighborhood  of  the 
singular  point  is  to  set  up  the  indicial  equation  and  find  its  roots. 
When  this  is  done  for  Equations  (3.2)  and  (3.3)  it  is  found  that  the 
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roots  are:  0,  0,  1/2,  1/2.  The  presence’  of  the  two • double  roots  indicate 
that  logarithmic  terms  must  be  included.  Thus,  each  variable  has  the 
form 


P1  +  (log  c)  P2  +  Q1/2  P3  +  U1/2'  log  O  P4  (4.1) 

where  the  P_^  denotes  power  series  with  unknown  coefficients.  Recurrence 
relations  obtained  by  substituting  these  forms  into  the  differential 
equations  determine  the  coefficients  of  the  power  series  in  terms  of  the 
arbitrary  constants. 

Since  a  numerical  solution  to  Equations  (3.2)  and  (3-3)  is  obtained, 
a  lengthy  discussion  of  the  singularity  is  not  needed.  Suffice  it  to 
say  that  the  constant  terms  in  the  series  P2  and  are  zero  for  each 
variable. -  Thus,  each  gradient  function  has  the  form 

1/2 

constant  .+•  (£  .  x  constant)  (4.2) 

in  the  neighborhood  of  C  =  0. 

5.  NUMERICAL  SOLUTION 

Equations  (2.12)  are  first  solved  separately;  then  Equations  (3.2) 
and  (3.3)  are  solved.  The  Runge-Kutta-Gill  method  is  used  to  integrate 
the  equations  numerically 

With  high  speed  computers  it  is  possible  to  study  the  solution 

empirically  by  carrying  out  a  "detailed  approach"  to  the  body;  that  is, 

to  systematically  decrease  the  interval  size  in  the  Runge-Kutta-Gill 

procedure  as  £  approaches  zero,  ending  the  calculation  at  an  extremely 

_ TO  ’  “  '  . 

small  but  finite  value  of  £(^10  ).  This  is  the  method  used  in  the 

’  *  ♦  ^  ,  ■  .■  ■  ^  p  ■  «. .  ’ 

present  work.  While  admittedly  not  elegant,  this  procedure  does  produce 
answers  in  reasonably  short  times,  and  a  posteriori  tests  indicate  that 
the  values  obtained  are  correct.  Thus,  unique  limits  for  the  gradient 
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functions  are  found  as  the  integration  interval  and  limiting  c,  approach 

1/2 

zero;  also,  the  calculated  gradient  functions  exhibit  the  J,  variation 
near  zero.  Equation  (4.2),  which  was  demonstrated  in  Section  4.  In 
addition,  the  calculated  solution  to  Equation  (3.2)  agrees  with  the 
known  value  (ft  =  4>/v)  obtained  from  Taylor-Maccoll  flow. 

6.  RESULTS 

In  the  case  of  cones  (where  =  0)  the  shock  curvature  reduces 

to  the  second  term  of  the  expression  in  Equation  (3.1l).  This  quantity 

was  calculated  for  a  wide  range  jof  Mach  numbers  and  cone  angles  for 

nitrogen  at  300°K.  Results  are  Ishown  in  Figure  3  (solid  curves).  Each 

curve  has  a  vertical  asymptote  occurring  at  a  cone  angle  for  which  the 

Taylor-Maccoil  flow  is  supersonic  at  the  shock  wave,  but  subsonic  at  the 

body;  K  changes  sign  here  and  then  decreases  in  absolute  value  with 
w 

increasing  0  .  The  qualitative  behavior  of  is  similar  to  that  of 
shock  curvature  at  a  wedge  tip,  shown  also  in  Figure  3*  Figure  4  shows 
the  variation  of  the  pressure  gradient  at  the  tip  of  a  cone;  the 
"Crocco  point"  asymptotes  appear  here  also. 

From  results  of  Millikan  and  White  (1963)  some  representative  values 
/  q^  \  the  normalizing  quantity  for  distances,  are:  '  q^  '  =  5.2, 

0.26,  and  0.034  cm  for  0^  =  30°, 1  40°,  and  50°,  respectively,  at 

=  10,  pm  '  =  1  atmosphere;  for  other  p^  the  given  '  q^  '  is  to  be 
divided  by  the  number  of  atmospheres. 

In  Figure  5  the  gradients  of  pressure,  temperature,  and  density  at 
the  tip  of  a  cone  are  given  as  functions  of  cone  angle  for  =  10, 

T  '  =  300°K  in  nitrogen.  The  density  gradient  is  unusual  in  that  it  is 
non-monotonic  and  changes  sign;  this  change  in  sign  takes  place  when 
the  flow  on  the  cone  surface  is  subsonic.  A  similar  type  variation  is 
found  for  density  gradient  at  the  tip  of  a  wedge  and  is  indicated  in 
this  figure;  here  the  change  in  sign  occurs  for  supersonic  flow  at  the 


wedge  surface.  Similar  behavior  for  the  density  gradient  in  plane  flow 
was  found  for  the  case  of  dissociating  air,  using  the  gas  model  and 
gradient  formulae  of  Spurk ,  et  al.  (1966). 

As  mentioned  in  Section  1,  the  present  work  has  application  in  the 
calculation  of  axi symmetric  nonequilibrium  flows.  In  previous  work, 
Sedney  and  Gerber  (1963),  an  initial  region  of  frozen  flow  was  assumed 
in  order  to  start  a  nonequilibrium  flow  computation  by  the  method  of 
characteristics.  This  lead  to  results  illustrated  in  Figure  6  for 
pressure  variation  on  the  body  in  a  typical  case  of  flow  over  a  cone. 

The  results  are  clearly  in  error  near  the  tip;  however,  a  curve  drawn 
through  the  points  further  downstream  can  be  faired  back  easily  to  the 
correct  pressure  at  the  tip,  as  evidenced  by  the  solid  curve.  This 
furnishes  an  indication  that  the  results  downstream  are  correct.  The 
gradient  calculations,  indicated  by  the  dashed  lines,  give  further 
assurance  of  the  validity  of  the  characteristic  computations  downstream 
by  demonstrating  that  the  solid  curves  can  he  extended  hack  to  the  tip 
pressure  with  the  correct  slope.  Gradient  values  can  also  be  useful  in 
initiating  nonequilibrium  characteristic  flow  computations. 

The  present  results  can  he  compared  with  those  obtained  by  an 
application  of  Dorodnitsyn ' s  integral- method.  South  and  Newman  (1965). 
This  latter  procedure  is  an  approximate  method  which  yields  algebraic 
equations  for  the  curvature  and  gradients  at  the  tip  of  a  cone.  Figure 
7  shows  the  comparison  for  two  Mach  numbers.  The  agreement,  in  general, 
is  good;  it  is  found  to  improve  as  the  angle  interval  6  -  0^  decreases. 

7.  CONCLUSION 

A  method  is  given  for  determining  the  shock  curvature  and  flow 
variable  gradients  at  the  tip  of  a  pointed  body  of  revolution.  The  use 

of  variables  introduced  by  Chester  simplifies  setting  up  the  calculations 
- 

A  characteristic  calaulation  of  one  example  with  an  initially  negative 
density  gradient  was  performed.  The  gradient  changed  sign  at  a 
distccnoe  of  3x  (initial  frozen  length )  and  thereafter  was  positive . 

All  other  flow  variables  were  monotonic. 
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•  CHARACT.  CALC. 


Figure  6.  Application  of  pressure  gradient  calculations  to  check  characteristic  calculations 


9 

i 


as  compared  to  use  of  polar  coordinates.  Although  the  specific  source 
of  departure  from  equilibrium  considered  here  is  vibrational  excitation, 
the  same  technique  can  be  applied  to  a  dissociating  gas,  with  the  same 
qualitative  results  expected. 

It  is  shown  that  the  shock  curvature  at  the  tip  of  a  pointed  body 
of  revolution  is  the  sum  of  two  terms:  i)  the  curvature  for  a  curved 
body  in  frozen  flow  and  ii)  the  curvature  for  nonequilibrium  flow  over 
a  cone.  Also,  the  Crocco  point ■ for  nonequilibrium  flow  occurs  at  the 
same  Mach  number  .and  cone  angle  combination  as  for  frozen  flow. 
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1  Professor'S.  I.  Pai 

Institute  of  Fluid  Dynamics  & 
Applied  .Mathematics 
University  of  Maryland 
College  Park,  Maryland  207^0 

1  Professor  R.  Probstein 

Massachusetts  Institute  of 
Technology 

77  Massachusetts  Avenue 
Cambridge,  Massachusetts 
02139 

1  Dr.  A.  Charters 

Defense  Research  Laboratories 
General  Motors  Corporation 
Santa  Barbara,  California 
93108 

1  Dr.  J.  Clark 

Division  of  Applied  Mathematics 
Brown  University 
Providence,  Rhode  Island 
02912 

1  Dr.  S.  Maslin 

Research  Institute  for  Advanced 
Studies 

The  Martin  Company 
7212  Bellona  Avenue 
Ruxton,  Maryland  21212 
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1  Dr.  A.  Puckett 

Systems  Development  Laboratory 
Hughes  Aircraft  Company 
Florence  Avenue  &  Teal  Street 
Culver  City,  California  90232 

1  Dr.  W.  Sears 

Graduate  School  of 

Aeronautical  Engineering 
Cornell  University 
Ithaca,  New  York  14850 

1  Dr.  R.  Sedney 

Research  Institute  for 
Advanced  Studies 
The  Martin  Company 
7212  Bellona  Avenue 
Ruxton,  Maryland  21212 

1  Dr.  J.  Sternberg 
The  Martin  Company 
Baltimore,  Maryland  21203 

1  Dr.  J.  Lukasiewicz 

Von  Harman  Gas  Dynamics 
Facility 

ARO,  Incorporated 
Arnold  AFS,  Tennessee  37389 

1  Mr.  H.  Nagamatsu 

General  Electric  Research 
Laboratory 
P.0.  Box  1088 

Schenectady,  New  York  13205 


1  Mr.  P.  Rose 

AVCO-Everett  Research  Laboratory 
2385  Revere  Beach  Parkway 
Everett,  Massachusetts  02149 

1  Mr.  J.  Whitfield 

Von  Karman  Gas  Dynamics  Facility 
ARO,  Incorporated 
Arnold  AFS,  Tennessee  37389 

1  Mr.  H.  Yoshihara 

General  Dynamics/CONVAIR 
P.0.  Box  1950 

San  Diego,  California  92112 

Aberdeen  Proving  Ground 

Ch,  Tech  Lib 
Air  Force  Ln  Ofc 
Marine  Corps  Ln  Ofc 
Navy  Ln  Ofc 
CDC  Ln  Ofc 
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